1. The spatial distribution of global Rs sites

Global spatial distribution of soil respiration sites

Global spatial distribution of soil respiration sites

Rs data from cold regions are more important, but how to increase the number of measurements? * Make equipment work normally in very cold and remote conditions * Increase resources devoted to Rs measurements * Methodological improvements: for instance, measuring once per day to get daily mean

Rs measured at diurnal soil temperature

Rs measured at diurnal soil temperature

Rs measured at mean annual soil temperature

Rs measured at mean annual soil temperature

2. The objects of this analysis are

3. Methods

Data

Statistics

Update Bahn’s model

3.2 test the relationship between Rs_annual and Rs_mat

## Fri Apr  5 16:11:37 2019  -------------------------------------------
## Fri Apr  5 16:11:37 2019  Bahn relationship for these data:
## 
## Call:
## lm(formula = Rs_annual ~ Rs_TAIR, data = sdata)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -745.35 -111.13  -39.72   86.33 1274.51 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 160.24086   13.26967   12.08   <2e-16 ***
## Rs_TAIR       0.94443    0.01625   58.13   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 194.9 on 821 degrees of freedom
## Multiple R-squared:  0.8045, Adjusted R-squared:  0.8043 
## F-statistic:  3379 on 1 and 821 DF,  p-value: < 2.2e-16

## [1] "first method-------------------------------------------test whether intercept=0 and slope=0 for [y-1x ~ x] model"
## [1] "p_intercept = 0, p_slope = 0.000657459332537873"
## [1] "second method-------------------------------------------test intercept and check whether get same results"
## [1] "t_inter = 12.076, p_inter = 5.07443172548757e-313"
## [1] "second method-------------------------------------------test whether intercept differ from 1"
## [1] "t_slope = 3.42, p_slope = 0.00065745933253803"

3.3 Ts sources (MGRsD, MGRsD_TAIR, From paper, Rs_Ts_relationship)

3.4 Annual Rs or Ts coverage effect

3.4.1 Test whether outliers affect the regression (need update)

  • Conclusion: need to update the code and results # moved to section 4.1

3.4.2 Effect of maximum allowed divergence between global climate data set and site-specific air temperature

  • Does TAIR_dev and TAIR_LT<_dev affect the relationship – YES!!!!!
  • TAIR_LTM_dev = with( srdb, abs( MAT_Del - MAT ) )
  • Does TAIR_LTM_dev () pull the slope off 1? – YES!!!!!
  • TAIR_dev <- with( srdb, abs( TAnnual_Del - Study_temp ) )
  • Figure E. Effect of maximum allowed divergence between global climate data set and site-specific air temperature, when given. As we throw out data points with high divergence, R2 goes up (top panel) and RSE goes down (bottom, g C m-2 yr-1).

3.4.3 Effect of maximum allowed divergence between annual precipitation from paper and Del

4. Results

4.1 Using Ts, TAnnual or MAT

4.1.1 Using soil temperature

4.1.2 Using T_Annual

4.1.3 Using MAT

4.2 Analysis when Rs_mat cannot represent Rs_annual

4.2.2 Does Ecosystem_type affect the relationship between Rs_annual and Rs_mat?

## geom_path: Each group consists of only one observation. Do you need to
## adjust the group aesthetic?
## geom_path: Each group consists of only one observation. Do you need to
## adjust the group aesthetic?
## geom_path: Each group consists of only one observation. Do you need to
## adjust the group aesthetic?
## geom_path: Each group consists of only one observation. Do you need to
## adjust the group aesthetic?
## geom_path: Each group consists of only one observation. Do you need to
## adjust the group aesthetic?
## geom_path: Each group consists of only one observation. Do you need to
## adjust the group aesthetic?

## geom_path: Each group consists of only one observation. Do you need to
## adjust the group aesthetic?
## geom_path: Each group consists of only one observation. Do you need to
## adjust the group aesthetic?
## geom_path: Each group consists of only one observation. Do you need to
## adjust the group aesthetic?
## geom_path: Each group consists of only one observation. Do you need to
## adjust the group aesthetic?
## geom_path: Each group consists of only one observation. Do you need to
## adjust the group aesthetic?
## geom_path: Each group consists of only one observation. Do you need to
## adjust the group aesthetic?

4.2.3 Does Meas_method affect the relationship?

4.2.4 RA- or RH-dominated effect?

4.2.5 Biome effect?

## geom_path: Each group consists of only one observation. Do you need to
## adjust the group aesthetic?
## geom_path: Each group consists of only one observation. Do you need to
## adjust the group aesthetic?
## geom_path: Each group consists of only one observation. Do you need to
## adjust the group aesthetic?
## geom_path: Each group consists of only one observation. Do you need to
## adjust the group aesthetic?
## geom_path: Each group consists of only one observation. Do you need to
## adjust the group aesthetic?
## geom_path: Each group consists of only one observation. Do you need to
## adjust the group aesthetic?

## geom_path: Each group consists of only one observation. Do you need to
## adjust the group aesthetic?
## geom_path: Each group consists of only one observation. Do you need to
## adjust the group aesthetic?
## geom_path: Each group consists of only one observation. Do you need to
## adjust the group aesthetic?
## geom_path: Each group consists of only one observation. Do you need to
## adjust the group aesthetic?
## geom_path: Each group consists of only one observation. Do you need to
## adjust the group aesthetic?
## geom_path: Each group consists of only one observation. Do you need to
## adjust the group aesthetic?

4.2.6 TAIR and precipitation variability effect?

4.2.7 Drought effect?

5. Discussion & questions

6. More analysis in the future